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Jacek Dziarmaga and Krzysztof Sacha 
Institute of Physics and Centre for Complex Systems, 
Jagiellonian University, Reymonta 4, 30-059 Krakow, Poland 
(Dated: December 08, 2005) 

In the Bogoliubov theory a condensate initially prepared in its ground state described by station- 
ary Bogoliubov vacuum and later perturbed by a time-dependent potential or interaction strength 
evolves into a time-dependent excited state which is dynamical Bogoliubov vacuum. The dynamical 
vacuum has a simple diagonal form in a time-dependent orthonormal basis of single particle modes. 
This diagonal representation leads to a gaussian probability distribution for possible outcomes of 
density measurements in position and momentum space. In these notes we also discuss relations 
with the U(l) symmetry breaking version of the Bogoliubov theory and give two equivalent gaussian 
integral representations of the dynamical vacuum state. 



INTRODUCTION 

It has been established in a number of papers 0, Q 
that a density measurement on a Bose-condensed state 
effectively "collapses" the state to a single Bose-Einstein 
condensate. The Bose-condensed state before the mea- 
surement is a quantum superposition over ideal conden- 
sates J T><fi il>((j)) \N : 4>) with different condensate wave 
functions 4>{x) but as position. (|14I15|1 s of more and more 
particles in the state are measured the state of remain- 
ing particles gradually collapses towards one of the ideal 
condensates \N : </>). As a result the density measure- 
ment outcome is p(x\4>) — N\(j>(x)\ 2 with probability 

p(0) « k/>(0)| 2 . 

In our recent paper Q we developed a measurement 
theory predicting probability of different density mea- 
surement outcomes on a Bose-Einstein condensate evolv- 
ing under external time-dependent perturbation. At zero 
temperature in the framework of the Bogoliubov the- 
ory the state of the condensate is described by a time- 
dependent dynamical Bogoliubov vacuum. The dynam- 
ical vacuum has simple diagonal form which leads to 
a simple gaussian probability distributions for different 
density measurement outcomes. The measurement the- 
ory for dynamical Bogoliubov vacuum is important be- 
cause in many experiments, like phase imprinting of dark 
solitons or condensate splitting in atom interferome- 
ters 0, manipulation of the condensate generates sub- 
stantial dynamical depletion which can qualitatively af- 
fect measured density patterns. 



DIAGONAL BOGOLIUBOV VACUUM 

According to Bogoliubov theory with well-defined 
number of atoms |g] a condensate which was initially 



prepared in its ground state (i.e. a stationary Bogoli- 
ubov vacuum) evolves into a time-dependent excited 
state which is formally a time-dependent or dynamical 
Bogoliubov vacuum annihilated by time-dependent Bo- 
goliubov quasiparticle annihilation operators. The dy- 
namical Bogoliubov vacuum can be brought to the diag- 
onal form 31 



f m \ 

44 + x <*°iai ) 

v a=l I 



N/2 



(1) 



with a time-dependent orthonormal basis of single par- 
ticle states (p a (t,x) and real time-dependent eigenvalues 
A a (t) € [0,1). In the diagonal state Q we sum over 
finite number M of non-condensate modes, keeping in 
mind that M — > oo. However, in practical calculations 
which almost always involve some numerics one is forced 
to work with a finite number of modes. This is why in the 
following we will keep finite M in all summations without 
any further comment. As the diagonal vacuum was de- 
rived elsewhere |3j , here we only mention that (f>o (t, x) is a 
condensate wave function solving time-dependent Gross- 
Pitaevskii equation and 4> a (t,x)'s are eigenmodes of the 
density matrix 

M 

m—1 
M 

^[l + dJV Q (f)] &(f, x)ct> a (t, y) . (2) 



Here [u m (t, x), v m (t, a?)]'s are Bogoliubov modes solving 
time-dependent Bogoliubov-de Gennes equations. In the 
number-conserving Bogoliubov theory modes u m (t,x) 
and (t, x) are orthogonal to the condensate wave func- 
tion <fio(t,x). dN a (t ) is average number of atoms depleted 
from the condensate wave function <f>o(t,x) to the mode 
4> a {t,x). The eigenvalues X a (t) in the diagonal vacuum 
GJ are 
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Phases of the modes cj> a (t, x) are chosen in such a way 
that the eigenvalues of the matrix 



A I 



E v * m { t ^) u m{t-,y) = 

m— 1 

M 

VdN a (t)[l + dN a {t)] <t> a (t,x)ct> a (t,y) ■ (4) 



are real and positive. For the sake of completeness we 
also quote here the diagonal density matrix for non- 
condensate modes 



A I 



Y v m (t,x)v* m (t,y) = 

m—l 

M 

^cUV a (t) 0*(t,2)0 a (t,0) 



(5) 



a=l 



One of the results of Ref . [3] is that the operators Ij2l4l5|l 
can be simultaneously diagonalized by the same single 
particle non-condensate modes <fi a (t,x). Phases of these 
modes are uniquely defined by the requirement that the 
eigenvalues of the operator (@J are real and positive. 



SUPERPOSITION OVER CONDENSATES 

Having in mind applications in quantum measurement 
theory it is useful to rewrite the diagonal state Q as a 
gaussian superposition over TV-particle condensates 

d qe -^*=i— q - \N : </>(x\q)) . (6) 



with the normalized condensate wave functions 



4>{x\q) 



M%) + -jjf E«=l qa<f>a(x) 



1 + 7? 2/3=1 ifj 



(7) 



Here the state \N : cf>(x\q)} is perfect condensate with all 
N particles in the same condensate wave function (f>(x\q). 
The (rather technical) proof of the equivalence between 
the diagonal vacuum in Eq.JQ) and the gaussian super- 
position in Eq.Q is given in Ref. 3]. The equivalence 
is approximate, it requires that the average number of 
depleted particles is much less than the total number of 
particles, 



A I 



dN = Y dN <* < N ■ 



(8) 



a=l 



This is not a new assumption but the usual requirement 
in the Bogoliubov theory. 



CONNECTION WITH THE SYMMETRY 
BREAKING APPROACH 

This section is a brief digression on the standard U{1) 
symmetry breaking version of the Bogoliubov theory 
where one splits the field operator into a c-number con- 
densate part plus a small quantum fluctuation 



#(x) 



(9) 



The quantum fluctuation is further Bogoliubov trans- 
formed as 



■ip{x) 



M 

E 

m—l 



(10) 



The initial stationary Bogoliubov vacuum state evolves 
into a time-dependent Bogoliubov vacuum state annihi- 
lated by the operators b m . In the U{\) symmetry break- 
ing Bogoliubov theory the reduced single particle density 
matrix in the time-dependent vacuum state is 



A I 



(x)tp(y)) = Y v m {t,x)i* m {t,y) 



and the anomalous density matrix is 

AI 

(iP(x)ip(y)) = Y u m (t,x)v* m (t,y) 

m—l 



(ii) 



(12) 



However, these density matrices are not confined to the 
subspace orthogonal to the condensate wave function 
(j>o(t,x) because, unlike in the number-conserving the- 
ory, here the Bogoliubov modes u m (t, x) and u^(t, x) are 
not exactly orthogonal to 4>a(t, x). This error can be cor- 
rected using the relations 



u m (x) = u m ($) - <p (x)((f> \u ri 
v m (x) = v m (x) - <t%(x)((j)o\v n 



Q% v m (x) (13) 



between Bogoliubov modes in the number-conserving and 
the symmetry-breaking theories After the projection 
on the subspace orthogonal to (j>o(t, x) we obtain the cor- 
rect non-condensate density matrix 

q% Q ff (ftffltm = 



E v m{t,x)v* m {t,y) = 

m=l 

oo 

J]d7V Q (i) ^* a {t,x)4> a (t,y) 

and the correct anomalous density matrix 

Qx Qg {^{x)Tp{y)) = 

OO 

^ u m (t,x)v^(t,y) - 



(14) 



m—l 

oo 

E 

a=l 



x)4>a(t,y) ■ (15) 
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The last equalities in Eas. (|14ll5[) follow from equations 
(|5I4|) in the number-conserving theory. 

Results of the conventional 17(1) symmetry break- 
ing theory can be translated into correct results of the 
number-conserving theory either by 

• projecting the Bogoliubov modes on the subspace 
orthogonal to the condensate wave function like in 
Eqs.CSJ), or 

• projecting directly the density matrices 
($(x)ip(y)) and (tp(x)4>( y)) in Eqs. flHT2l) 
in a way described by Eas. (|14ll5() . 

or any other way to elliminate contamination by the con- 
densate wave function cf>o(t,x) of what is wrongly be- 
lieved to be purely non-condensate density matrices. 



BOGOLIUBOV REPRESENTATION 

We have shown that in the limit of large N the two 
forms Q and JSJl of the Bogoliubov vacuum are equiva- 
lent. However, there is yet another representation of the 
vacuum which may appeal more to some of the readers 
because it is very similar to the coherent state represen- 
tation of the vacuum in quantum optics. The represen- 
tation is 



|0i 



d A b e 



\N : <}>(x\b)) 



(16) 



Here d 2 b = UZ^i d 2 b m , b*b = £^ =1 b* m b m and <j>($\b) is 
a normalized condensate wave function 



<j>(x>\b) = 



(17) 



1 M 

■y= [ b m U m {t,x) + b* m V*Jt,x) 

v m— 1 



where the normalization factor is 

A I 



N 1 = l + ^ E {b m u m + b* m v* m \b n u n + b* n v* n ) . (18) 



m.n—l 



In this Section we will prove equivalence between Eq. 1(16(1 
and Eq.© in the limit of large N. The main difference 
between these equations is that in Eq. (JSJ there are M in- 
tegrations over real coordinates q a while in Eq. i|16f) there 
are also M integrations but over complex coordinates b m . 
In this sense Eq.@is more compact representation than 
Eq.(j]3J|. 

As a first step we use the large N limit to transfer the 
normalization factor Af to the exponent in Ea. l|16fl . 



d 2 b 



M 

E 

TCI — 1 



[K 



|0) . (19) 



Here the new operators are u m — (u m \ip) and 
v m = {v* n \ij)). At this point we make a 
transformation to the orthonormal basis of 4>aS, 
Em=i [ b mvL + b^v^ =] J2a=i or equivalently 

Y \bmUrn + K n V* m ] = ^ Za<f>a • (20) 
m a 

After this transformation the state 1|19|) becomes 



i M \ N 

— ^z.ai] |0). (21) 

v a—l / 



Here 6's are linear functions of z's 

M 



b m = ^ [ Z "( U ™l^a) ~ Z a{ V m\<t>* a )} ( 2 2) 



a=l 



obtained by inverting the transformation l(20|) . Using this 
inverse transformation and equations IL'1 1151) we rewrite 
the exponent in Eq.JSlJ as 

J2 b * m b m +J2z* a z a = 

rn a 

5^2(i + diV >*« Q - 



y/dNjl + dNa)z*zl - y/dN a {l + dN a )z a z a (23) 

With this exponent the gaussian integral in Eq. 1)21(1 gives 
a simple correlator for z's 



d 2 z e ^£ m h '" fcm ^£= 2 ° ; 



1 + diV, 



P 



(24) 



In the last equality we use the relation in Eq.©. The 
gaussian state in Eq.J2U is completely determined by 
this correlator. The same correlator for real g's is also 
obtained in the real gaussian integral 



dq e 



-*£„ 



9/3 <Z 7 = Spy Xp 



(25) 



Comparing the correlators 124fl and (|25|l we find that the 
gaussian state (|21|) is equal to the state 

r 2 ( i M \ N 

e-lE.lt U + ^E^J 1°) ( 26 ) 

with real coordinates q. Finally after normalization of 
the creation operator we get 



dq e 



E, 



I % + ^- E a =l 9a«L 



A? 



i + n E/3=i ip 



|0)(27) 
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and this state is the gaussian representation of Bogoli- 
ubov vacuum in Eq.©. In this way we have shown that 
the representations © and (ft!))) are equivalent. 

The complex representation with Bogoliubov coef- 
ficients b m has twice as many integration variables as the 
real representation © with real coordinates q a . What 
is more, the real representation © is spanned by the or- 
thonormal basis of </> a 's while the complex representation 
is spanned by the Bogoliubov modes (u m ,v m ) which are 
not orthogonal in the usual sense. The real representa- 
tion reveals the diagonal structure of the vacuum which 
is implicit in the complex representation. 



PROBABILITY DISTRIBUTION 

Having shown equivalence between different represen- 
tations of the dynamical Bogoliubov vacuum (|1I6I16|I 
we now return to the main subject of this paper which 
is probability distribution for different condensate wave 
functions. The quantum superposition in the real repre- 
sentation © directly leads to the probability distribution 
for different q's 



P(q) cx [] e 



(28) 



which is valid for the non-condensate modes with large 
occupation numbers d7V Q when we can replace 



1- Aq 

2X n 



1 



4dW, 



(29) 



compare Eq.@. The distribution P(q) is the desired 
probability distribution for different possible density 
measurement outcomes given by 



p(x\q) = N\^(x\<l)\ 2 



N 



<fo(*! x) + -j= E a =i q»0 a (t, X) 



1 



N 2^/3=1 1/3 



(30) 



Equations (|28|l and (|3U[1 define a simple scheme how to 
simulate possible density measurement outcomes on the 
dynamical Bogoliubov vacuum state. 



MEASUREMENTS IN MOMENTUM SPACE 

Up to this point we assumed ideal density measure- 
ment where all particle positions are measured at the 



same time. However, there is no reason why these results 
should be limited to measurements of particle momenta 
instead of position. In fact the diagonal second quantized 
representation prefers neither position nor momen- 
tum representation, and consequently all the following 
formulas can be rewritten in momentum representation 
simply by replacing x — > k and y — > p. In the end we 
arrive at the probability distribution P(q) in Eq.JSHJ) for 
different possible outcomes of momentum distribution 



p(k\q)=N 



(j>o{t,k) 



Y,a=l1a<t> a (t,k) 



IV" a 2 

N 2^/3=1 1/3 



(31) 



A quantum measurement on the dynamical Bogoliubov 
vacuum state which measures all particle momenta at the 
same time will give momentum density distributions i|31|) 
with probabilities l|28|) . 



CONCLUSION 

In conclusion, we derived a convenient diagonal form of 
the time-dependent Bogoliubov vacuum which greatly fa- 
cilitates simulations of quantum measurements on Bosc- 
condensed systems both in position and momentum 
space. 
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